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Abstract 



This paper presents a new approach to studying galactic structures. They are considered as the low- 
^ ' frequency normal modes in a disc of orbits precessing at different angular speeds. Such a concept is an 

adequate alternative to the commonly used approach of treating the disc as a set of individual stars rotating 
at near-circular orbits around the centre. The problem of determining the normal modes is reduced to a 
CO ' simple integral equation in the form of the classical eigen-value problem, where the eigen-value is directly 

equal to the pattern speed of the mode, Qp. An examination of the general properties of the basic integral 
equation shows that two types of solutions exist, bar-like and spiral. The numerical solutions of both types 
are obtained. The characteristic pattern speeds are of the order of the mean orbit precession speed, although 
for the bar-modes flp can markedly exceed the maximum precessing speed of orbits. It is shown that the 
' bar-mode grows due to the immediate action of its gravitational field on the stars at the resonance regions. 

1^-^ , As for the spiral mode, its excitation is probably due to the inner Lindblad resonance that can promote 

' mode growth. 



Cf^ , Keywords: Galaxies: structure. 

o 

'p^! 1 Introduction 

O I A disc galaxy is primarily a set of stars which rotate around the centre at near-circular orbits at the angular 

■4^ ■ velocity r2(r). The observed spiral and bar- like structures are customarily treated as perturbations in such a 

^ I differentially rotating disc. Thus, one might expect that the typical pattern speeds of these structures, Op, 
should be of the order of some average star angular velocity, Q. In reality, however, the pattern speeds Op are 

. ^ ^ of the order of the characteristic precession speed of star orbits, Opr, which is only a small fraction of O. Recall 

i that Opr = 0(r) — K{r)/2 for near-circular orbits, where «:(r) = (40^ -I- dSl^/dr)^/^ is the epicyclic frequency. 
H ' According to Lynden-Bell (1979), if stars involved in the formation of the structures satisfy the inequality 

e= |0p-0p,.|/0< 1, (1) 

then each star orbit as a whole, but not the individual stars, participates in the perturbations^. Consequently, 
for studying such structures, it is more reasonable to use the model of the disc of star orbits precessing with 
different speeds than the commonly used concept of the differentially rotating disc of individual stars. The main 
conjecture of the present paper is that the galactic spirals and bars are the normal modes in such a model of the 
disc galaxy. In view of its simplicity, this approach, as we will see below, allows one to clarify the underlying 
physical mechanisms in the formation of galactic structures. 

The Lynden-Bell inequality mentioned above are justified by results of numerous calculations (e.g., 
Athanassoula & Sellwood (1986) for bar-modes or Lin, Yuan & Shu (1969) for spirals). 

Fig.QJi shows 0(r) and Opr(r) for the Plummer potential '^o{r) = — (1 + r^)~^/^ used by Athanassoula & 
Sellwood (1986) in their A^-body bar-mode analysis. The solid thin horizontal line O™™ = 0.14, the dash-dotted 
line Op = 0.21, and the dashed line fl'p^^ = 0.3 correspond, respectively, to the minimum, mean, and maximum 
pattern speeds from the list given by Athanassoula & Sellwood (1986) in their Table 1. Fig.QJ) represents the 
ratios 50/0 ((50 = |0p — Op^l) for these modes. For the first mode localized inside the circle r = 2 (thin dashed 
lines in Fig. ^), typical ratios are of the order of 0.1. For the second mode (mean pattern speed), these ratios 

* E-mail: epolyach@inasan.rssi.ru 

^Note, after Arnold (1989), that it was Gauss who had proposed, for studying the perturbations of planets by each other, to 
smear out a mass of each planet along its orbit in proportion to a time and replace the attraction of planets by the attraction of 
such rings. 




Figure 1: Justification of the Lynden-Bell model of precessing orbits for bar-modes studied by Athanassoula & 
Sellwood (1986): a) the curves Q{r) and flpr{r) for the Plummer potential (thick solid lines), and the minimum, 
mean, and maximum patterns speeds (thin solid, dash-dotted, and dashed lines, respectively) from the list of 
Athanassoula & Sellwood (see their Table 1); b) the ratios SU/^l for fip from a). 
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Figure 2: The same as in Fig.^ for the model of the Galaxy from Lin, Yuan, & Shu (1969). 



are typically of the order of 0.2. Even for the fastest mode {flp = typically, Sft/ft ^ 0.3, taking into 

account that this mode is more concentrated to the centre. Therefore it is little wonder that the pattern speeds 
calculated below (Section 5) in the framework of our approach repeat very accurately those of Athanassoula & 
Sellwood (1986) in all cases, including the most rapid bar-modes. 

On the other hand, it is evident from the curves in Fig. ^ that the inequality ^ and hence our approach 
fail at sufSciently large radii. This is certainly of no importance for determining the pattern speeds, since the 
bar- modes are mainly localized near the centre. Roughly speaking, the ratio (5f2/ri becomes of the order of 1 
near the corotation radius (i.e., at Vc ~ 3.6, 2.7, and 2.0, for the first, second, and third modes, respectively). 

Fig.|3L,b show similar graphs for the classical spiral example from Lin, Yuan & Shu (1969) where the 
model of our Galaxy was considered. The pattern speed used here corresponds to the mean of the interval 
(rip = 11 — 13km s^^ kpc^^) suggested in the cited paper, i.e., the horizontal line in Fig.|2t is flp = 12km s^^ 
kpc^^. Fig. shows that the analysis of the eigen- modes of Galactic disc for this model could be performed in 
the framework of our approach (without the restrictions of the WKB J theory) . Note that substantially higher 
values of fip, compared to those from Lin et al. (see, e.g., Blitz 1983), are most likely due to an independent 
peripheral stage of the spiral structure. Such structures can be analyzed using the integral equations described 
in the end of Section 2. 

We picture the disc as comprising a large number of slowly deforming and rotating elliptical rings, along 
which an individual stars moves so fast it is no more than a blur. Small changes with radius in the shape or 
orientation of rings give rise to regions of enhanced surface density (see, e.g., Fig 6-11 of Binney & Tremaine 
1987). In the absence of self gravity any initial pattern of over-densities will wind up because the part of the 
pattern at radius r will pecess at a rate flpr that varies significantly with r. 

Rings that have inclined minor axes exert gravitational torques on one another. These torques cause the 
angular momenta and radial actions of the rings slowly to change, thus deforming the rings and altering any 
initial pattern of surface density. The torques also affect the precession rates of the ring's the minor axes. The 



interaction between rings is strongest when the rings are physically adjacent, i.e., have similar values of the 
actions. If we neglect the non-negligible range of the interaction, the disc has something in common with a goods 
train of deformable wagons connected by springs. Self gravity allows waves of ring deformation to propagate 
through the disc just as the springs permit waves of wagon displacement and deformation to propagate down 
the train. We show that these waves can generate a pattern that rotates rigidly at an angular frequency that is 
larger than the fastest precession rate in the disc. 

The picture has points in common with Lynden-Bell's (1979) theory of bar formation in that it focuses on 
slowly precessing orbital rings, but it does not conclude, with Lynden-Bell, that the bar rotates at a rate that 
compromises between the fastest and slowest precession rates of the consituent rings. Instead we find that both 
bars and spirals, in common with many wave phenomena, move faster than the underlying medium. Several 
studies have shown that the rotation rates of bars in N-body models conform to our predictions rather than 
those of Lynden-Bell. 

The model of precessing orbits can be introduced directly as is done in Section 2 after Eq. ((T7|) . However, 
this model is rigorously justified by the use of the regular procedure of the perturbation theory in the small 
Lyndcn-BcU parameter e from (1), with averaging over fast radial oscillations of stars. This programme is 
implemented in the beginning of Section 2. The resulting integral equations for disc normal modes derived in 
two ways, namely, (i) directly in the model of precessing orbits, and (ii) with the help of averaging, are identical. 

Lynden-Bell (1979) showed that the quantity J/ = Ji -I-/2/2 (/i and I2 are the usual actions) is an adiabatic 
invariant, provided that the inequality ^ is satisfied, i.e., for sufficiently slow perturbations. This means that 
the orbits evolve by the action of the potential of low-frequency modes only through change in the angular 
momentum L, with Jf held constant. 

Assume that in the initial moment the distribution function of orbits is To{Jf,L). Let AL be a small 
variation of the orbit angular momentum owing to the action of a small potential perturbation for the time from 
to to the current moment t. The flow of the phase fluid is incompressible, so the phase element with the initial 
density jFo( J/, L — AL) should come into the point (J/, L) at the moment t. Accordingly, the Euler differential 
of the distribution function, i.e., its perturbation, is equal to 

T = MJf.L - AL) - MJf.L) ~ -J-^AL (2) 

(JFq = dJ^al Jf, L)/dL). T depends on the derivative T'q. So it is not suprising that just this derivative is found 
to be critical to our theory of disc slow modes. Henceforth, this derivative is referred to as the Lynden-Bell 
derivative of the distribution function^. It turns out that the behaviour of the Lynden-Bell derivative determines 
the type of eigen- modes (bar-like or spiral). 

To develop the above considerations, immediately resulting from the remarkable paper by Lynden-Bell 
(1979), into the theory of low-frequency disc modes, one should invoke the dynamical equations and the ex- 
pression for the 2D disc potential (see Sections 2,3). As a result, we obtain the integral equations for these 
eigen- modes (see Eqs. H12|l and (jlfifl below). If a potential satisfles these integral equations, the orbits evolve 
in such a way that the pattern rotates uniformly, with a certain angular speed fJp. The physics of our bar- like 
and spiral solutions consists just in the fact that these are (slow) normal modes. All specificity (compared to 
general modes) is concerned with their slowness. 

The ensuing sections of this paper present the following material. In Section 2, we derive the basic integral 
equation of the theory (by two methods). The basic integral equation is then reduced to the form of the classical 
eigen- value problem. Section 3 is devoted to the analysis of the general properties of the basic integral equation. 
In particular, the crucial role of the Lynden-Bell derivative of the distribution function is revealed. In Section 4, 
this derivative is examined in more detail by the example of the Schwarzschild model. In two following sections, 
the specific solutions of the basic integral equation are given. To demonstrate the capabilities of the proposed 
theory, test models studied earlier by the A^-body method (Athanassoula & Sellwood 1986) are adopted (Section 
5). The theory gives results in close agreement with results of the A^-body simulations. A wide variety of spiral 
modes are studied in Section 6. In the last Section 7, a brief summary of the results is given. 

2 Basic equations 

The most convenient variables that should be used under studying the low- frequency modes are J/ = /i -I-/2/2 
and L = where Jf is the Lynden-Bell adiabatic invariant (Lynden-Bell 1979), while (/i, 12) = I are the usual 
actions. The angle variables corresponding to I are w = (wi,W2). The variables {Jf,L) = J are also actions, 
for which w = (u;i,iD2) = {wi,W2 — wi/2) are the canonically conjugate angles. One can show that for the 

^When studying the slow bars, Lynden-Bell (1979) has used the same derivative of the star precession speed, dflpr/dL. 
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npri3)t, (3) 

where t is the time, rtpr{J) = ^l2(J) — ^^i(J)/2 is the precession speed of the orbit with the actions J, fii and 
^2 are the frequencies of star radial and azimuthal oscillations: fli — dH{){i) / dli [i — 1,2), = v^/2 + $o('') 
is the unperturbed star energy expressed in terms of I, v is the star velocity, $o ('') is the equilibrium potential. 
In the limit of near-circular orbits Qi = n{r), U,2 — ^{r)^ where k is the epicyclic frequency. 

Eq. means that the angular variable W2 is slow because \^pr\ <SC |^^2|- It is easy to show that iv2 — a, 
where a is the azimuth of the orbit's minor axis (see below, Fig.|2I). 

For later use we give here the coordinates of stars r, W2 — ^p, expressed in terms of J, wi. The radius 
r = r{J,wi) is determined by solving the equation 

r 

r dr' 

w,ir,3)=n, / , (4) 

(it is assumed that < wi < tt). The slow angular variable 

W2 ^ (p + (piiJ,wi), (5) 

where 

r 

dr' 



2E{J) - 2$o(r') - 



The distribution function /(J,w, i) is governed by the collisionless Boltzmann equation 

% = [HJ], (7) 

where H is the Haniiltonian of a star in the self-consistent gravitational field (/)(r,t), 

1 



[H, f] denotes the Poisson bracket. 



\H fl = , _ dj^dH_ _ dj^dH 

^ ' ^ ~ dwi dJf dw2 dL dwi dJf dw2 dL' ^ ' 

To linearize the collisionless Boltzmann equation, we assume that 

/ = .Fo(J)+.F(J,u)i)e^("'^^-"*), 
= $o(r) -I- $(J, wi)e*("'^^-'^*\ 

where .7-o(J) = /o(I) and JF(J, wi)e'*^"""2~'^*^ are the unperturbed and perturbed distribution functions, re- 
spectively, $(J, wi)e*'-"™^~"*'' is the perturbation of the potential, m is the azimuthal wavcnumbcr, lu is the 
frequency (a; = mQp). As a result, we obtain the equation 

^t{u;-mnpr)T + ni—^—^—+im^—^. 9 
owi oJf owi oL 

To find the desired low- frequency solutions, one can use the perturbation theory in the small Lynden-Bell 
parameter e from Let T = JF'^' + JF'^) + . . . be the perturbation series in powers of e, so that T^^^ is 
obtained from © by neglecting the terms proportional to fip^ and <I> oc G: dJ^^^^ / dwi = 0, i.e. J^^^^ = J-"'^-'(J) 
is an arbitrary function of the integrals of motion, which is subsequently specified by using the periodicity 
condition for the solution of the next approximation^. 

''Since \lj — mQpr \ <C f2i, the wave frequency oj is not determined in the first approximation. 



A 



The equation for JF^^^ takes the form 



-i{u;-mnpr)J'^^> = —^— + im<i>—f. 10 

owi oJf owi oL 

Given the periodicity of functions JF^^-' and <&, averaging over wi in the interval (0, 2tt) yields 

271- 

-(cc;~mr!p,)^(i) «m^$ ($ = ^ ^ dwi $). (11) 



The azimuthal number m can take only even values. Formally, it follows from the relation exjp(imw2) = 
exp(imw2)-exp(— ir7ii(;i/2), which remains periodic in wi (with a period equal to 27r) only for even m. Physically, 
a closed precessing orbit is an oval symmetric relative to the centre. So the torque from the perturbations with 
odd m will break rather than rotate such orbits'*. 

Calculating the perturbed surface density, 

and using the expression for the 2D disc potential, we find 



$(r) ^-G j dr'^^ = 



ri2 

where ri2 = [r^ + r'^ — 2rr' cos{lp' — Lp)]^/"^ . Now one can change in the last formula from r', v' to J', w' taking 
into account that dr'dv' = dJ'dw' (dJ' = dJydL', dw' = dw']^ dwj ) : 



$(J,wi) ^-G j dJ'dw 



, _ , JF(i)(J') exp[imi5w2] 



7-12 

with 5w2 = w'2 — W2- Finally, averaging the potential $ over iZ;i,we obtain the following integral equation^; 

$(J)^|/dJmJ')^^;;^^W, (12) 
where T'^ = {dWf, L)/dL)\j^ = dh/dh - l/2dfo/dh, 

n(J, J') ^ / d«)idu)ld5«)2^^^P^^^. (13) 
J ri2 

The function 11 can be reduced to the following suitable form: 

TV TT 

n(J,J')=8 J dOui cosirvpi J dw[ cos ■mip['ip{r,r'), (14) 



where 



, , , , cos ma , , 

^(r,r') = / da (15) 

\/ r-^ + r' — 2rr' cos a 







The function n(J, J') (multiplied by the factor imG/2n) has a meaning of torque acting from the harmonical 
distribution of orbits, e*™"'2^ with the action J' to the orbit with the action J and orientation ■W2 — 0. 



■^It is obvious that the self-consistent perturbations made by such ovals must repeat treir symmetry ip ^ ip + it, which is valid 
only for modes with even m. Among them, the bi-symmetric mode m = 2 is predominant. For instance, just this bar-mode is most 
likely the only unstable mode. Indeed, the excees of Op over (Op7.)max is due to the self-gravitation. This effect is proportional to 
the function 11, as it is seen from I16i . From I14i . it follows that the function 11 decreases with m. So, the growth of the modes 
with m > 4 owing to CR and OLR would be much loss than for the m = 2 mode. The ILR resonance for m > 4, however, will 
prevent the growth. 

^The integration over real J' is correct only for unstable frequencies (Imuj > 0). Thus roots with (Imo; < 0) should be omitted. 



It is notable that the integral equation (|12|) in terms of the function JF^^) has the form of the classical 
eigen- value problem, where the eigen- value is directly the pattern speed flp. Indeed, determining $ through 
!F^^^ from Hll|l and substituting it into (|12|l . one obtains 

npT'^^\3) = J d3'K{3,3')T(^\3'), (16) 

with the kernel 

K{3,3') = ^^^(J)n(J, J') + npr{3)6[3 - J']. (17) 

Eq. H16|) is the basic integral equation of the theory. 

The integral equation H16|l is much simpler than the general integral equations for disc normal modes derived 
earlier by Kalnajs (1965) and Shu (1970). An ordinary desktop PC is quite capable of solving the integral 
equation lfTC|l . 

The most important advantage of the present approach is the fact that it makes clear the underlying physical 
mechanisms of the instability processes developing in the disc. To reveal these physical mechanisms by using the 
integral equations of Kalnajs or Shu would be very difficult, and the same is true for the TV-body simulations. 

The formal derivation of the integral equations (|12II16|I for the low-frequency modes of gravitating discs 
given above is rigorous. However, the physical situation will be more transparent after deriving the equivalent 
integral equation in another way using explicitly the fact that each orbit as a unit should take part in slow 
perturbations of interest. 

Accordingly, let us consider the distribution function of closed precessing orbits, f{Jf,L;a,t), such that 
dJA — fdJfdLda is the mass of stars at orbits within a given interval dJfdLda, a is the azimuth of minor 
axis, so that the orbit precession speed Qpj.{ Jf, L) = a. The coUisionless kinetic equation for such a distribution 
function is 

dt ~ dt + + 9L 



where one takes into account that {df/dJf)Jf — for the slow modes of interest (Lynden-Bell 1979), and 
L = M, where M is the torque acting on the orbit with given Jf,L,a. Note that such a distribution function 
/ and kinetic equation (|18|l were earlier suggested by V. Polyachenko (1992). 

By linearizing the equation 118(1 and assuming that /i, $i cx exp(— zcji + ima), one can obtain the equation 
that is analogous to Ea. lfTT|l : 

-i{u: - m%r)h = (19) 

where Mi is the perturbation of the torque. Let us check that the equations H19|) and (|11|) are actually identical. 
First of all, it is easy to show that 

M, . (20) 

where $i is the potential averaged over the selected orbit; this potential is produced by all other orbits of a 
system: 



I 



1 



$1 = $i(J/,L;a) = f- / dsci>i(r)p|:;f^^)(s), (21) 



where pun{s) = l/v{s) is the linear mass density characteristic for each orbit {v is the total star velocity, s is 
the current length of the orbit, r is the radius-vector at the orbit), fi(Jf,L) = /dsp/i„(s) is the mass of stars 



on the orbit, 

f ds'da'dJ'dL'p,:^J'^'\s')fi{J'L'; a') 
$i(r) = -G/ / I'hn V JJn _l ^22) 



{ip and ip' are the current azimuths at two orbits). 

In the formulae 1(20(1 and ((21(1 . we used the same notation $ for the averaged potential as earlier in the 
derivation of the integral equation 1(12(1 . since in both cases it was actually the same quantity. To convince 
oneself that this is correct, one can immediately compare the two expressions for taking into account that a 
star is within the interval ds during dt — ds/vtot = dr/vr — dwi/fli. Thereafter it remains to make sure that 



the slow angular variable W2 and the azimuth a of the minor axis are identical. Indeed, jSJl can be rewritten in 
more informative manner: 

W2=^ (fi- Aip + flprAt, (23) 

where At is the time it takes for the rotation of a star through the angle Aip between the azimuth of the 
minor axis and the current azimuth (p. Then the identity ^2 = a is clear from Fig.O After invoking Poisson's 
equation, we obtain the integral equation that is coincident with Eq. ifT^ . 

In addition to the derivations above, let us indicate the simplest way to use Eq. Q of the Introduction for 
obtaining the relation Ijlll) . The canonical Hamilton equations in variables (J,w) are 



dH ■ dH 
dH , dH 



where H{Jf, L;iui,W2) is the Hamiltonian. Since the angular variable W2 is "slow", for studying the low- 
frequency modes, the procedure of averaging the equations of motion over a quick variable (in this case, Wi) is 
appropriate. As a result, we obtain (see, e.g., Arnold 1989) 

27r 

1 f OH 
Jf ~ 7^ I dwi— = 0; 



27r J dwi 



I f dH a$i 



271 J ^ dw2 dw2 ^ 



The first equations means the adiabatic invariance of J/, while the second equation determines the evolution of 
the angular momentmn L. Assuming that $ corresponds to a normal mode, from (|25|1 one can find: 

Integrating the last equation over t from to = ^00 to t, we find 

AL = g-»(c.-mOp.)t^ (26) 

taking into account that the perturbation is switched off at t — > — cx). Substituting H2t)|l into Q we obtain 
It is clear that the perturbed precession speed is the same periodic function of time as AL. 

Recall that the integral equation H16() is valid only for even m. However, one can obtain the integral equation 
describing the low- frequency modes for m divisible by mo (where too is an arbitrary integer) in the form (jl6|l . 
when replacing 

Wl o ^1 dfo 1 dfo 

W2 = W2--, npr^n2--, — = — --— 

with 

W2 ~ UJ2 , i 'pr — ^'2 , — or or' 

Too Too uL 012 riiQ oii 

respectively. 

The case too = 1 is particularly interesting. The corresponding basic integral equation is suitable for studying 
low-frequency modes of discs around black holes, in planetary rings, and so forth. Since only two types of the 
potentials, (x and oc are available, in which all orbits have zero precession speeds, there are only two 
cases (toq — 1 and toq = 2) where our approximation that wi is fast is valid in an extensive region. 



3 General analysis of the basic integral equation 

As it follows from 1)14(1 . the function H is real and symmetric: n(J, J')* = H(J, J'), H(J, J') = n(J', J). Let us 
restrict ourselves to the regular solutions for the averaged potential $. Taking into account Hill) , one can safely 
divide both sides of ((16|l by JFq, multiply by J-'^^^ , and integrate over J. Calculating the imaginary part of the 
resulting equation, one finds 

(ImOp). j dJ^^ = 0. (27) 



Figure 3: Illustration of the relation H23|l . which proves the identity W2 = a. Dotted and solid ovals show 
orientations of an orbit at times t = and t = At, respectively; the orbit is turned by an angle QprAt. The 
azimuth of a star initially (at t — 0) coincides with the azimuth of the minor axis (point A) . The position of the 

star at t = At is shown by point B; its azumuth is ip =xB (dash-dotted arc), and the variation of the azimuth 

is A(p —AB (dashed arc). 



By using Ea. (|27|l can be reduced to the form 

(Imrjp)-L™ = 0, (28) 

where denotes the angular momentum of the mode that can be obtained from the general expression of 
Lynden-Bell & Kalnajs (1972) retaining the only term which dominates for the low-frequency modes: 

r ij^'i'p r il'p 

'--l'"-Tr-J''^^w=k?- '''' 

It is found that the type of solution of the basic integral equation depends crucially on the behaviour of the 
derivative J-'q. Below, we consider two cases. 

1. Assume that J-q is strictly positive everywhere in the phase space of the system. Then Lm is negative, 
as well as the energy of the mode SEm, since Em = ^pLm (Lynden-Bell & Kalnajs, 1972). Hence, from (|28|> . 
one obtains Imilp — 0. The corresponding real eigen-functions, J-^^\ describe non-spiral solutions. We identify 
them with the inner parts of bar-modes, i.e., bars themselves. 

Obviously, in the case J^q > under consideration, the integral equation (|16|l determines merely the pattern 
speed Re flp of a mode. The question now arises of whether there are some external reasons for growth of bar- 
modes. It is found that the growth can be due to the exchange by the angular momentum with the resonance 
stars at the corotation and the outer Lindblad resonance (OLR)^. These resonances are located rather far from 
the central regions where the modes are mainly localized. In this respect these resonances are external relative 
to the modes. Of course, in the end the correctness of our approach is justified by the comparison with results 
of A^-body simulations (see Section 5). 

The corresponding growth rate can be calculated by the formula 

7 = (30) 

where Lm = Lm"^ + L^m , and the expressions for the rates of exchange by the angular momentum at CR (im'' ) 
and OLR {L\^>) can be obtained from the general formulas of Lynden-BcU & Kalnajs (1972) through their 
minor transformation: 



^If the disc is immersed into the real (but not rigid as in the paper by Athanassoula & Sellwood (1986)) spherical component, 
one should take into account, generally speaking, the resonance angular momentum exchange between the bar-mode and stars 
of these component. The dynamical friction caused by resonance interactions of stars of spherical systems with a wave was first 
studied by Polyachenko & Shukhman (1982). 



where S is the Dirac delta- function, ri(')(J) = fl2(J) + ^ l)f^i(J)/2, and 



= dJ'dwidu}'iV(^^')-^^'nJ')e~''"'^''^^'^"^'^ (32) 

is the Fourier coefHcient corresponding to CR {I ~ 1) and OLR (/ = 2) in the expansion of the potential in the 
series in e*'™^ . 

The interaction of stars in the resonance regions with the gravitational potential of the mode leads to the 
spiral responses. For the corotation resonance, as an example, the response consists of two parts. The resonance 
response of the disc, without taking its self-gravity into account, produces spirals of angular length about 7r/2 
(see, e.g., Polyachenko 2002a). The self-gravity also extends the resonance spirals by 7r/2 (Polyachenko 2002b). 
Therefore, the total angular length is about tt, which is typical for the majority of SB galaxies (see, e.g., Sandage 
1961). 

It should be emphasized that above we consider only the linear stage of the formation of structures in 
SB galaxies. The observed galaxies can have strongly nonlinear bars. The length of such bars is apparently 
determined by 4:1 resonance. This is exactly the case, e.g., for the nonlinear bar- modes computed by Sellwood 
& Athanassoula (1987). For a symmetric bar which is free of odd Fourier harmonics, in particular m — 3 
harmonic, the resonance 4:1 is the nearest to the centre, and thus is the most important one. 

2. Let us say that J-q becomes negative in some regions of the phase space. For the realistic distribution 
function (see Section 4), these regions, if they exist, can occupy only a small fraction of the total volume of the 
phase space. Then, in addition to the fast bars described above, the spiral-like solutions can occur. Contrary to 
the bars, these new modes grow due to their intrinsic instability, which is inherent to the mode itself. As follows 
from H28(l . for an unstable mode (Im flp > 0), the angular momentum — 0, i.e. the contributions to Lm from 
the regions with opposite signs of the Lynden-Bell derivative !Fq cancel each other exactly: = + = 0. 
Obviously, the instability is due to the inner Lindblad resonance, since the criterion of instability, JFq < 0, 
coincides exactly with the Lynden-Bell & Kalnajs condition of wave exitation at ILR (see (jHlll when I ~ 0). 
This resonance becomes the source of spiral waves in the case under consideration, contrary to the commonly 
accepted view of the role of ILR. Accordingly, the growth rate of the unstable mode can be calculated as 

^=2LT = ^' (33) 

where 

i± = ^ y dj.F^|$|25[f7p,(j) - rip], (34) 

r± 

r+ and r_ denote the regions of the phase space with positive and negative values of the Lynden-Bell derivative 
J-Q, respectively. Note that below the mode growth rates are obtained directly from the solution of the integral 
equation Hlt)|) as 7 — 2ImJ7p. However, it is also useful to derive a simple formula for estimate of the growth 
rate from (|33l I34() . The simplest case is that with one narrow domain of the negative derivative J-q. Using the 
latter equality of (|33|l and approximating the integrals as the products of the average values of the integrands 
and small volume of the phase space domain, we obtain: 

7-|f^;.|AL, (35) 

where AL is the width of the region with JFq < 0, and fi^^ = dilpr/dL is calculated at the narrow region of 
the phase space of interest. This estimate is consistent with the explanation of the physical mechanism of the 
angular momentum exchange at the ILR given by Lynden-Bell & Kalnajs (1972). However, two more points 
need to be made: (i) the sign of the effect is opposite to that in the cited paper, in accordance with the fact 
that J-Q < 0, and (ii) strictly speaking, their explanation is appropriate only for the case of near-circular star 
orbits, although, as we demonstrate below, the very elongated orbits play an essential role. Nevertheless, their 
considerations can readily be generalized to such a case. 

4 Examination of the Lynden-Bell derivative by the example of the 
typical model 

Let us consider the generalized Schwarzschild distribution function (Shu 1970) 



where E is the star energy, vq is the radius of the guiding centre: L = rQi}{ro), Ec{ro) = i;g(ro)/2 + <I>o('"o) 
is the star energy at the circular orbit, voiro) — roO(ro) is the circular velocity, ^o{ro) is the equilibrium 
potential. The specific model is given by the functions cro(ro) and co(ro); in the epicyclic limit, when vq/cq 3> 1, 
o'o{ro) — J^oiro), co{ro) = Cr{ro), where So(^o) and Cr{ro) are the surface density and radial velocity dispersion 
respectively. In the general case, Eo(''o) and Cr{ro) are expressed by aoiro) and co(ro) in a more complicated 
manner. 

The Lynden-Bell derivative of the distribution function (|36ll is equal to 

~ :^o<ra—, — -r-rQ— — - + Vq- 



dL rlK?{ro) \ f^(ro) K(ro) cro(ro) 

+ , J^., , + r.'-^iE - E^ir,))] , (37) 
co(ro) 2moS2(ro)ci^(ro) cg(ro) J 

where = /o(i?(J), ro(L)), a prime denotes the derivative with respect to rg. As a rule, J-"o > either in 

all the phase space or at least for the most part (almost everywhere). It is provided by the term r^K'^ /2raoVlc^ 
in 13 7() . For example, in the case of the flat rotation curve (fo = const), this term equals (\/2/'7io)(wo/co)^- 
Since uq/co 3> 1 almost everywhere in disc galaxies, this term is dominant. 

However, there may exist some narrow regions, in which such predominance is broken. Firstly, this is 
possible in the regions where the function ^(r) is sufRciently small. For instance, with a rotation curve similar 
to that of our Galaxy, such regions can be located at r « 2.5 kpc and 8 — 10 kpc. Note the obvious causes of 
such peculiarities in the rotation curves: (i) a changeover from the potential of the spherical component to the 
potential of the disc, and (ii) a sharp edge of one of the disc components. Secondly, the function co(ro) increases 
rapidly when approaching the galactic centre, so that the rotation velocity uq in the central parts becomes of 
the order of cq. 

Relatively small deviations in the behaviour of the rotation curve may lead to very different solutions of 
the basic equation. Some other factors are also very important. The excess of the pattern speed above the 
maximum precession rate is due to self-gravity. Thus the solutions have the bar-like form only in sufRciently 
massive discs, while for discs with relatively small masses we can obtain only spiral modes. The number of very 
elongated orbits plays the essential role. For example, the excess of such orbits over the number prescribed 
by the generalized Schwarzschild distribution function H36|) strongly increases the growth rate of the spiral 
modes due to increase of Ifip^l fo^' these orbits. Besides, it is seen from (|37ll that |JFq| (JFq < 0) increases with 
{E~E,{ro)). 



5 Numerical test solutions of the basic integral equation 

The basic integral equation is studied numerically. The unknown function !F^^\J) as well as the kernel K{3, J') 
are considered on the 31 x 31 grid in the phase space [E, L). Then the obtained matrix equation can be solved 
using standard methods of linear algebra. 

It is advisable to put aside the distribution function (|36|) until the next section and demonstrate the capa- 
bilities of the theory by applying it to test models. For such a role, we choose two models investigated by the 
A''-body simulations of Athanassoula & Sellwood (1986). The equilibrium potential for all the models is the 
Plummer potential, $o(^) — —(1 + r'^)~^^'^- 

First of all, we are dealing with a Kalnajs (1976) model, which Athanassoula & Sellwood (1986) have denoted 
by (to = 6, /3 = 0, g = 1, Jc = 0.25) (see their Table 1). The computed spectrum of the matrix equation, which 
consists of 31^ eigen- values, is shown in Fig.0] The frequencies are real, since the Lynden-Bell derivative for the 
Kalnajs models is positive everywhere in the phase space. Just few of eigen-values correspond to the discrete 
spectrum, while others mimic the continuous part of the spectrum^, which is the van Kampen - Case modes 
(van Kampen 1955; Case 1960). They are located in the interval between the minimum and maximum values 
of the precession speed: (—0.051,0.126). These modes are the stable solutions of the form oc (5(J — Jo). The 
peculiar van Kampen modes are of no interest in studying the bar modes. Thus, in the iV-body simulations 
of Athanassoula & Sellwood, the bar is always one of the modes of discrete spectra. The latter are the eigen 
solutions of the basic integral equation. Both the perturbed potential and surface density of the discrete modes 
differ in the number of radial nodes, the nodeless mode being the mode with the maximum pattern speed. 

In the spectrum given in Fig.0| one can see 5 discrete modes. The figure shows that the pattern speeds of 
the first and second modes (counting from the right) coincide with the pattern speeds obtained in the A^-body 

'^The existence of the continuous part of the spectrum is a property of the kernel I17i . which defines a bounded noncompact 
integral operator (c/. Reed & Simon 1972). 
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Figure 4: The spectrum of pattern speeds (filled circles) computed as the eigen- values of the problem for 
one of the Kalnajs models. The crossess are the complex pattern speeds, with the growth rates estimated by the 
formulas H3U|I . The diamonds are the "experimental" complex pattern speeds according to Athanassoula 

& Sellwood (1986). The arrows point to the minimum and maximum values of the star precession speed. 
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Figure 5: The spectrum of complex pattern speeds (filled circles) computed as the eigen- values of the problem 
((TCll . for the model {m = 6, (3 = 3, q = 1, Jc = 0.6) of Athanassoula & Sellwood (1986). The cross is the 
complex pattern speed, with the growth rate estimated by (|31ll . The diamonds are the "experimental" complex 
pattern speeds according to Athanassoula & Sellwood. The arrows point to the minimum and maximum values 
of the star precession speed. 

simulations to within 6%. It is found (see below) that the growth rates of the other 3 modes are much lower, 
so they could not be observed in the numerical experiments. 

Rough estimates of growth rates for the most rapid modes, in which we assumed that the star orbits are 
circular, give 71 — 0.117, 72 = 0.054. Each of these values is composed of two parts corresponding to corotation 
(CR) and OLR: jicr = 0.0199, jiolr = 0.0968, 72Cfl = 0.0107, 720Li?. = 0.0437. Note that in this case, the 
OLR terms are dominating. The mode with the maximum pattern speed has the highest growth rate. Obviously, 
the reason is that such a mode possesses the smallest corotation radius, where the disc surface density is high. 

The estimate for the growth rate of the most rapid mode is nearly twise as much as one obtained by 
Athanassoula & Sellwood (1986) in the iV-body simulations. For the second mode, the "experimental" and 
"theoretical" growth rates are quite close. Among the reasons of the possible discrepancy, one can mention the 
use of the circular orbit approximation. Recall that the estimates obtained by Athanassoula & Sellwood (1986) 
within the swing amplification approach are twice as little than their "experimantal" values. 

Another model was denoted by Athanassoula & Sellwood as (m = 6, /3 = 3, g = 1, Jc = 0.6)^. The 
computed spectrum of eigen-values of bi-symmetric modes for this model, is shown in Fig.El The continuous 
spectrum of van Kampen modes is again located in the interval between the minimum and maximum values of 
the precession speed: (—0.096,0.13). 

Since the Lynden-Bell derivative for the model under consideration is not positive everywhere in the phase 
space, the discrete spectrum contains complex modes. In Fig.jSJ we can see one unstable complex and one real 
discrete frequencies. The complex eigen-frequency obtained from the solution of the basic integral equation is 
very close to the "experimental" frequency, given by Athanassoula & Sellwood (1986). The spiral pattern of 
the unstable mode is given in Fig.lHI 

The pattern speed of the real mode coincides with the pattern speed obtained in the A^-body simulations 
to within 3%. The estimate of its growth rate, in the circular orbit approximation, gives 7 = 0.018, the OLR 

**The pattern speeds for other models from the Ust of Athanassoula & Sellwood (1986) were also computed as the eigen-values 
of the problem 1161 (about ten models in all). In all cases the results coincide very closely (with an accuracy of 10% or better). 
Besides, the eigen-frequencies of the basic integral equation for the isochrone potential coincide within a few percents with those 
obtained by Kalnajs (1978) using his matrix equation. 



Figure 6: The spiral pattern of the unstable mode in the model {rn = 6, P = 3, q = 1, Jc = 0.6) of Athanassoula 
& Sellwood (1986). The dashed circle indicates the position of the inner ILR (the outer ILR, for this model, is 
at r = 5). 

contribution being dominating. This estimate is 30% larger than that obtained by Athanassoula & Sellwood 
(1986) in the A^-body simulations. 

The models explored above are rather artificial. Nevertheless, they were convenient for us because it was 
thoroughly studied in the A^-body simulations by Athanassoula & Sellwood (1986). Besides, the analysis above 
points to the mechanism of bar-formation in the general case. It is clear that bars can form in galaxies with 
sufficiently massive discs. In this case, the necessary inequality Qp > (ripr)max can be fulfilled, otherwise the 
inner Lindblad resonances occur. Then, in the case of a positive Lynden-Bell derivative, only the continuous 
spectrum of the van Kampen waves exists. Note that the fulfillment of the inequality fip > (f^pr)max is made 
difficult not only for low-mass discs but also for discs with a high concentration of mass at the centre. 

The permitted frequency (pattern speed) of the galactic bar must be one of the eigen-values of the basic 
integral equation. It means that the corresponding waves remain unchanged over many revolutions of the galaxy, 
despite the different rates of orbit precession. The half-turn spirals adjacent to the bar form due to the same 
resonance interaction. 

6 Spiral solutions of the basic integral equation 

Here we consider in greater detail the situations that arise when the Lynden-Bell derivative is negative in some 
regions of the phase space. Then a great number of available possibilities appears; they are briefly outlined 
below. An extended discussion will be published elsewhere. 

As a base distribution function for the analysis of the spiral solutions of the basic integral equation, we choose 
the generalized Schwarzschild distribution function described in Section 4. For the standard model, we assume 
the "Milky Way-type" rotation curve (see Fig.[7^), and the exponents cro(^o) — aQe~'^°/^'', Co(ro) = CQe'^"^"^" that 
fix the particular Schwarzschild function. In our units, one kpc is equal to 1, = 3, = 2rii, the gravitational 
constant G = 1, (Tq = i'^fd)^^- The magnitude of the rotation curve is chosen to match the disc component 
of the velocity curve at radii vq > 5. The maximum rotation velocity is (wo)max = 1-2, cq ~ 0.83, so that for 
rg = 8 the ratio i'o(''o)/co(''o) — 4. For such parameters, the galactic disc is exponential almost everywhere, 

r 

excluding the very central part. Given the rotation curve, the equilibrium potential is $o('') = J dr' VQ{r')/r' . 
In such a manner we obtain the model of a galaxy, which is similar to our Galaxy (but, of course, the model 
galaxy is not identical).^ 

Fig-Et shows the Lynden-Bell derivative, J-q, for the standard model calculated by (|37|) at the line of circular 
orbits {E = Ecr{ro)). There is only one narrow region with JFq < located near the centre. Fig.|Sl demonstrates 

^Remind that the occurrence of the grand design spiral structure in the Galaxy is so far open to question. The galaxies, which 
we mean in our theory, are more in the nature of NGC 2997 (see, e.g., the image of this galaxy on the cover of the book by Binney 
& Tremaine (1987). 




Figure 7: (a) The "Milky Way- type" rotation curve (solid line); dashed line corresponds to the contribution of 
the exponential disc; (b) The precession speeds, fipr, of the circular orbits for the standard model (solid line) 
and the model with a "peak" near r = 8 — 10 (dashed line). The horizontal straight lines correspond to different 
pattern speeds; (c) The Lynden-Bell derivatives, !Fq^™'°\ for the circular orbits. 



how deep into the phase space the negative values of JFq penetrate. 

The basic integral equation for the standard model provides a large number of unstable modes. They can 
be either trailing or leading spirals, the latter are the majority. Narrowly-localized van Kampen solutions of 
the continuous spectrum, which occupy the wide band of real frequencies, fi™" < < ^^pr^, can be both 
the trailing and leading spirals since for each of these modes, effectively, the gravitational constant G 0^°. 
When employing the concept of eigen-modes, just the packets of leading van Kampen waves participate in 
the classical swing amplification process described in the original papers by Goldreich & Lynden-Bell (1965), 
Julian & Toomre (1966), and Toomre (1981). If the dispersion of orbit speeds, AQpr, is sufficiently small (that 
is always assumed), we obtain quasi-monochromatic packets of van Kampen waves. Just such wave packets 
moving with some group velocity from the central regions and approaching the corotation, can be amplified 
by the resonance effects under the transformation from the leading form to trailing one. Though each of the 
van Kampen modes is essentially kinematic (G — > 0), the finite packets of these waves are governed by their 
self-gravity. In particular, the local dispersion relation of Lin-Shu-Kalnajs should be valid for them (away from 
the resonances). 

Here we touch on those improvements, which the present theory can introduce into the original scheme of the 
swing amplifier. They are connected with the possibility of unstable discrete leading spiral modes. In particular, 
there is then no need to reflect the short trailing waves somewhere in the central regions to close the feedback 
loop. Given that the discrete leading modes are monochromatic, the efficiency of the resonance amplification 
should grow in comparison with quasi-monochromatic packets of the van Kampen waves. These modes might be 
the "embryo" waves for their subsequent resonance growth due to the swing amplifier. Then, apparently, from 
the full set of these waves, the mode with the maximum potential perturbation at the corotation is selected. 
Let us emphasize again the "inversed" role of the ILR (compared to the commonly accepted one), which acts 
as the generator of spiral waves. 

The growth rates of spiral modes for the standard model are rather low: they provide the growth of the 
amplitude e times, typically for (3 — 5) • 10^ years (taking into account that (i'o)max = 1-2 corresponds to 
the rotation period T = 2.5 • 10^ years). The main reason for such low values of growth rates follows from 
Fig.lHl This figure shows that the curve Qpr{E,L) = Reflp (for a typical unstable spiral mode) intersects 
the strip of negative values of the Lynden-Bell derivative, !Fq, deep inside the phase space of the system 
under consideration, at {E, L) corresponding to highly elongated orbits. But the number of such orbits in the 
Schwarzschild distribution function (decreasing exponentially with the ratio {E — Ec{ro))/cQ{ro)) is small. The 
growth rates of the spiral modes can be significantly increased if one takes into account that the number of such 
orbits is actually very much larger. Indeed, we are dealing with the central regions of spiral galaxies, where the 
contribution of the spherical component, i.e. stars of Population II that have the radially elongated velocity 

^"As computations show, the van Kampen modes are the sole type of solutions for the case when > everywhere in the 
galactic disc with the assumed parameters. 
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Figure 8: The strip of negative values of the Lynden-Bell derivative, T'q, at the phase plane {E,L). The phase 
domain is bounded by the line of circular orbits (from the top) and the line corresponding to 0.02 • (Jx))max, 
where (J-o)max is the maximum value of the distribution function, J-'q{E, L), at each fixed value of E (from the 
right). The solid curves, npr{E,L) = fip, indicate the locations of resonance orbits on the phase plane, for 
several values of Or,. 



diagram, is essential. It implies that the spherical component can play an important active role, instead of the 
passive role of a rigid halo as is commonly accepted. Computation with the standard model that is accordingly 
modified, leads to tenfold growth rates. The modification consists in the substitution CQ^ro) — > 2co(ro) provided 
that {E — Ec{ro)) > c§(ro). For such a modification, the observed surface density and radial velocity dispersion 
are almost unchanged. It is also important that this modification leads to the predominance of trailing spiral 
modes. The observed galactic spirals are most likely just these modes (or similar spirals in other modifications 
of the standard model). For instance, the modes are most often the trailing half-turn spirals as in the majority 
of observed galaxies. 

In particular, an interesting modification consists of forming the "peak"^^ on the curve ^Ipr at sufficiently 
large r, where the Lynden-Bell derivative can be either positive or negative (in Fig.|Sl the upper right corner of 
the {E,L) phase plane). We omit descriptions of the variety of cases, restricting ourselves to the only picture 
of the typical spiral mode (Fig.EJ. 

In conclusion let us touch on one very important aspect of the problem. Fig.[3: shows that the "generalized" 
Lynden-Bell derivative J^q^"^°^ at mg > 2 can be made negative more easily compared to the case mg = 2. This 
effect is due to the inverse dependence of the large positive term r^K^ /2moilcQ in the expression 1)^^7(1 on niQ. 
Hence, in the outer parts of galaxies (for our Galaxy, just beyond the Solar circle) the additional multi-arm 
structures are generated, which are independent from the central spiral structure (the latter usually has two 
arms)^^. Furthermore, one can expect that several modes will be generated there simultaneously leading to 
formation of the so-called flocculent structures (Elmegreen & Elmegreen 1982, 1987). All the aforementioned 
agrees with observations of galaxies (see, for example. The Hubble Atlas of Galaxies, Sandage (1961)). As a 
rule, in the central part of the disc, the regular two-arm spiral is seen, which is then superseded by a multi-arm 
or disordered structure. It is appropriate to mention here that Elmegreen (1992) has clearly demonstrated that 
three-arm spiral structures appear only on the periphery of galaxies beyond the inner Lindblad resonance for 
the m = 3 mode. 



7 Discussion 

The new approach to the study of structures in spiral galaxies suggested above is based on a view of a galactic 
disc as a set of precessing orbits. Following this model, we have obtained the simple integral equation that is a 
very convenient tool for studying galactic structure. Non-axisymmetric structures are considered as the normal 

^'^For our Galaxy, such a "peak" has been marked since the well-known paper by Lin, Yuan, & Shu (1969). 

^^The curves J^q^"^"'^ (r) in Fig.|7|; are shown only at sufficiently large radii to emphasize that modes with mg > 2 cannot typically 
occupy as wide a region as the two-arm mode does. This is explained by the different behaviour of flpr{r) and f2p™'''(r) for mo > 2. 




-10 -5 5 10 



Figure 9: One of typical pictures of unstable spiral modes. 

modes of the basic integral equation H16|l . It is found that the properties of solutions depend crucially on the 
behaviour of the Lynden-Bell derivative, !Fq, at the phase space of the disc. The solutions have bar- like or spiral 
form depending on whether J-q is positive everywhere or negative somewhere in the phase space. 

Within the framework of the suggested approach, which does not include the swing amplifier, we succeeded 
in calculating the pattern speeds and growth rates of all the modes obtained earlier in the -/V-body simulations 
by Athanassoula & Sellwood (1986), the coincidence of results being very high. 

Thus, it is shown that the bar-mode grows due to the immediate action of its gravity on the stars in the 
vicinity of corotation leading to the angular momentum exchange between these stars and the bar-mode. Note 
that Athanassoula & Sellwood (1986) could readily calculate exact growth rates of all the modes they computed 
by the formulas (|30l I31() . However, they preferred the language of swing amplification for rough estimates of 
the growth rates. 

To avoid confusion, we should emphasize that here we deal only with the standard fast bars. As for the 
slow bars of Lynden-Bell (1979), they are most likely the central parts of the unstable spiral modes (in this 
sense, they are secondary). Note that here we leave aside the possibility of the slow bar formation in very hot 
centres of galactic discs (in full analogy to the ellipsoidal deformation of spherical systems under the action of 
the radial orbit instability). In contrast, the fast bars are primary while the adjacent spirals are secondary since 
the spirals form as the response of the disc near the corotation to the bar gravity. 

The commonly used approach to the problem of formation of the normal spiral structures consists in using 
either the swing amplification mechanism (Toomre 1981) or its weaker counterpart (overreflection with reference 
to the waser mechanism - see, e.g., Bertin & Lin 1996). Note that the swing amplification mechanism is 
sometimes the only way for growth of the initial perturbations. It seems likely that such galaxies do not have 
any organized spiral structure of the mode nature. Certainly, this mechanism can effectively amplify transient 
wave perturbations. 

In traditional mechanisms, the region of corotation plays the central role. In our approach, we draw attention 
to ILRs that can under certain conditions cause the excitation of a great variety of spiral modes. Unlike the 
usual epicyclic idea of ILR as a certain circle, ILR in our unstable modes extends over wide regions of the 
galactic disc. This is due to participation of a variety of orbits including highly elongated orbits (see Fig.lS)). 
On the one hand, ILR can give rise to leading spiral modes, which in turn can be the "embryos" for further 
work of the swing amplifier. On the other hand, ILR can provide a variety of trailing spirals including those in 
the region between two ILRs (see Fig.[7|D). 
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Appendix. Some additional arguments in support 
of the used approach 

1. The perturbed distribution function / is represented as a Fourier series that contains resonance denominators: 

^ (uj — mih) + /sii 

It is natural to name the terms with / — —1, I — 0, and ^ = 1 as "ILR", "CR", and "OLR" terms, respectively. 
In Section 2, we have restricted ourselves to the only ILR term. As it was shown in Section 4, this leads to the 
correct results (with an accuracy of 10%) in determination of the pattern speeds. Thus, one can conclude that 
the ILR-term has a dominant role. 



2. The ratio of the denominator of the ILR term to the denominator of the CR term (closest to ILR) is 

\Lj-mn2\ \n-np\' ^ ' 

where the quantity SQ — — mflpr\ was already used in the Introduction. Strictly speaking, this ratio is equal 
to the Lynden-Bell small parameter e from ^ only when \n\ 3> \^p\, i.e. in the central regions far enough from 
the corotation. Fig.llUI shows the ratio (|39|l for the same 3 models as in Fig.^ K is seen that A ^ 1 for the 
slowest mode, so the validity of the theory is justified for such modes. The fact that the theory has a certain 
field of application justifies the main conclusions of the theory (classification of modes, physics of instabilities), 
derived in Section 3. They can be adequate even beyond the framework of rigorous validity of our approximate 
theory. For the other two modes from Fig.^J the ratio is considerably smaller than 1 only in the central regions. 
Nevertheless, we can hope that the approach is still valid even for these modes, at least quantitatively, taking 
into account that the perturbation amplitude decrease with radius, being relatively small at the corotation. 
Besides, there is one more factor that weakens the role of CR term; it is considered in the next point. 



3. The coefficient ai in the expansion (|38(l is proportional to the linear combination of the derivatives of the 
unperturbed distribution function /o(I) (see below, item 5): 



where 



a/o(i) , 9/o(i) 



dh dh ■ 

Thus, for the bi-symmetric mode m — 2 one should actually calculate the ratio 

/^,02(I) ,/^,-12(I) 



A' 



(40) 
(41) 

(42) 



pr\ 



rather than H39|) . In the epicyclic approximation, which is valid everywhere beyond the very central parts of 
galactic discs: 

(43) 

Fortunately, for the CR term I — 0, so the ratio (|42|l contains the additional small factor (compared to H39|l 'l. 
Apparently, this may be very important for explaining why our theory is valid even for fastest modes. 



dfo 




dfo 


dh 


> 


dh 



4. Intuition suggests that higher harmonics (with \l\ > 2) cannot considerably influence the large scale modes, 
which are the subject of our study. The particular number of harmonics, needed for the mode computation, is 
determined "experimentally" . The comparison with A'^-body simulations shows that the only ILR term provides 
a reasonable accuracy. In the next point we give a more detailed analysis of the relative roles of the terms in 
the expansion H38I) . 
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Figure 10: The ratio A{r) for different values of the pattern speed Vlp-. sohd, dash-dotted, and dotted Unes 
correspond to the minimum (f^™™ = 0.14), mean {Clp — 0.21), and maximum (fi™'''^ — 0.3) pattern speeds, 
respectively. 



5. Let us derive the generalization of the basic integral equation that takes into account all a resonanses of 
the perturbation with fixed azimuthal number to. Let the perturbed potential and the perturbed distribution 
function are expanded in Fourier series 

^{\w,) ^Y.^i{l)e''^\ /(I,u;i) = ^/z(I)e'''"\ (44) 

where I = (/i,/2) are usual radial and azimuthal actions, Wi is the radial angular variable. Substituting the 
full potential and the full distribution function 

cf>o + <i>(I,^«l)e^'"'"^ /o + /(I,ii;l)e""'"^ (45) 

into the collisionless Boltzmann equation, one can obtain after linearization 

/i(/r!i+TOl72-c^) = a>i/^,,„, (46) 

where the function /g ^.^^^^ is 

klmW^l^J^+^^J^- (47) 

The second relation between (I) and /; (I) is provided by the Poisson equation 

/f(r' \pimw'^ 
drdw^ ^^ ' ^' . (48) 

Multiplying both sides of the equation by e""™!/^ and introducing function 

ipi{wi) = W2- wi/2 - if (49) 

(see Section 2), one can obtain 

^$KI)e'^'+'"/^^'"' j dI'du;i^/r(I')e'(''+™/^^'"^e™'''^iV(^,^')- 



(50) 



I 

where Sipi — ip\ — ipi. Denoting 11; ;/ the integral 

n;,;,(I,I') = y" dwidw'iV(^?'')e'^''^™^^^'"^"*^'^™^^^""e™'^'^\ (51) 
one obtains the second relation: 

<i>z = -^ 101^^,(1,1')/;' (I') (52) 

(here and below the summation over the blind index V is assumed). Thus, expressing /; from (|46|l one can 
derive the final set of integral equations: 

$Ki) = -^\ drnu,(i,iO/o,r,n(iO ,o — (53) 



For numerical calculations, it is more convinient to rewrite the equation H53|) in the form of the classical 
eigen-value problem. Indeed, using H46I) one can obtain 



Mi){m, + mn,-u;) - -^/^,,„(i) J di'n;,,,(i,i')//'(i') (54) 

Using the discretization of the integrals J dl' ^ AI' , the eigen- frequencies can be obtained from the condition: 

G 



det 



2^ mw (I, I') AI' + E(Zf7i (I) + mn, (I) - u;) 



0, (55) 



where E is the unity matrix. Note that the basic integral equation H16() is the special case of H54II when 
1 = 1' = -1. 

In fact, the equations H54II55|I give a new formulation of the general eigen-value problem for galactic discs. 
This formulation can be considered as an alternative for the well-known matrix approach of Kalnajs (ApJ, 
V.212, p.637 (1977)). 

Let us employ the equation H55|l to define more exactly the eigen-frequencies of the bi-symmetric m = 2 
modes in the Kalnajs model (6, 0, 1.0, 0.25), which was already studied in Section 4. Athanassoula & Sellwood 
reported that there are two unstable bar-modes: uJi = 0.465-|-0.066z, uj2 = 0.33 -I- 0.058i. We found that for this 
model the Lynden-Bell derivative is positive everywhere, so only bar-modes with Op > (ripr)max were possible. 
Using the basic integral equation IjlGII . we can determine only the pattern speed of the modes; they found to 
be Rewi = 0.44, Rew2 = 0.35. We argued that the growth rate for these modes can be obtained, taking into 
account the angular momentum exchange at the resonances (primarily, CR and OLR). Rough estimates gave 
us 71 = 0.117 and 72 = 0.054; notably, the OLR contribution were dominating in both cases. 

Restricting ourselves to consideration of three terms with \l\ < 1 in the set of equations (jHU, we can study 
the impact of three main resonances: ILR, CR, and OLR. First of all, we obtained the unstable bar-modes, 
when all three terms are present (run 1). Then, the influence of a particular resonance can be "measured" 
separately by omitting an unnecessary term. In such a way, we study the influence of CR term, switching off 
the OLR term (run 2). After that we study the impact of the OLR term, omitting CR term (run 3). In all these 
three cases, when the principal ILR term is present, the pictures of the unstable modes in the complex w-plane 
resemble each other. The situation has changed when we dropped the ILR term. In the result we obtained 
unrealistic picture of a sea of unstable modes with growth rates less than 7 • 10~^. This fact certainly confirms 
the dominating role of the ILR term. 

The computed eigen-frequencies are given in the table. 



Run 


Resonances 


Mode 1 


Mode 2 


1 


ILR, CR, OLR 


0.48 + 0.058i 


0.38 -1- 0.024i 


2 


ILR, CR 


0.43 + 0.015i 




3 


ILR, OLR 


0.49 -1- 0.036i 


0.38 + 0.024i 



It can be seen that: 

• All pattern speeds gained small positive biases, compared to the pattern speeds obtained from the basic 
integral equation H16|l . 

• When all three terms are taken into account, both the pattern speed and the growth rate of the most 
unstable mode (mode 1) is determined with very good accuracy. Separate consideration of CR and OLR 
(runs 2 and 3) shows that the contribution of OLR term into the growth rate exceeds one of CR. 

• From comparing the eigcn-frequency of mode 2 obtained in runs 1 and 3, it follows that mode's growth is 
entirely due to the interaction of the gravitational potential of the mode with the resonance stars at OLR. 

In principle, equations (|54|l allow to consider any number of terms, although it can be checked that the 
terms with |Z| > 2 are of less importance. For example, taking account of terms I = ±2 will not lead to any new 
unstable modes. The fastest modes become rather more unstable: 71 = 0.08, 72 = 0.032. 
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